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In this article I summarize results on global observables, hadron spectra, and ratios of 
integrated hadron yields as presented at the Quark Matter 2002 Conference. I also attempt 
to put these results into context and convince the reader that an evolving coherent picture 
begins to form, shedding light on the state of matter created in relativistic heavy ions 
collisions. However, at this conference we have been presented with such a wealth of new 
qj | data on hadronic signals that no summary can give proper credit to everyone. A definitive 
summary of all the results on hadron spectra, ratios, and yields and their interpretation 
would certainly require a much longer paper. I therefore limit myself to a compilation of 
results that I consider to be the most interesting, and refer the reader to the large number 
of excellent papers given at the conference for further details. 

It was obvious that this conference was dominated by the latest findings from the RHIC 
experiments. Despite the short time span between the last run and the Conference, spectra 
of many particle species were reported, both over a wide range in and, for the first time, 
also systematically over a broad range of rapidity. The longer RHIC run in 2001/2002, 
delivering for the first time the design energy of y/s NN = 200 GeV, allowed the experiments 
to accumulate more statistics and, together with the progress in the understanding of the 
detectors, resulted in higher precision measurements than possible at the lower energy 
(\/ s nn = 130 GeV) studied in the previous run in 2000/2001. 

Although the SPS heavy ion program at CERN has been almost completed, with all 
experiments delivering the physics they were built for, work on the analysis of data, their 
systematic comparison and interpretation continues. In order to understand the physics 
behind the state of matter produced in relativistic heavy ion collisions one must, more 
than everything else, understand the excitation functions of all observables involved. Only 
if we are able to describe the evolution and behavior of the systems as we increase the 
energy - from AGS, through SPS to RHIC - will we gain insight into the rich variety of 
physics we are facing. 

In this article I start with a summary of global observables presented at this confer- 
ence followed by a discussion of particle ratios and chemical freeze-out conditions. I 
then address the question of boost-invariance at RHIC and conclude with a discussion of 
transverse momentum spectra and kinetic freeze-out parameters. 
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One of the earliest probes suggested for QGP formation involves a study of the global 
parameters of the events, e.g. the energy deposition, multiplicity, and the average trans- 
verse momentum of the emitted particles, as a function of center-of-mass energy y/s NN , 
mass number A, and centrality of the collision. For example, by studying the multiplicity 
of the produced particles one might estimate theoretically the entropy produced in the 
collision. Sudden changes in behavior with varying centrality or A would be indicators of 
a phase transition. So far, however, no such anomalous changes have been observed, at 
either the AGS, SPS, or RHIC. All results on global observables shown at this conference 
indicate a rather smooth evolution in centrality and ^s NN . This, of course, does not 
necessarily imply the absence of a phase transition, but might be rather an indication 
of either the insensitivity of these observables to the early phase of the collision and/or 
might suggest a second order phase transition (or a cross-over). 

With the commencement of the RHIC 
program the question of multiparticle pro- 
duction in nuclear collisions became more 
complex due to the poorly understood role 
of perturbative QCD (hard processes). 

The study of charged particle multiplic- 
ity as a function of the number of par- 
ticipating nucleons was conducted by all 
four RHIC experiments and the results pre- 
sented in Fig. |l] depict the high quality 
and agreement between the experiments 
|]. It should be noted that, even more 
than in the determination of the multiplic- 
ity and the required corrections for decays, 
absorption, and feed-down, the difficulties 
of this analysis lie in the determination of 
Npart, i-e. the extraction of the underly- 
ing collision geometry from the data. This 
requires a detailed understanding of the 
trigger, especially trigger efficiencies, con- 
tamination and the study of possible auto- 
correlations. All four RHIC experiments 
now use Monte-Carlo Glauber model cal- 
culations similar to the ones implemented 
in the Hijing model as compared to nu- 
merical calculations of the nuclear overlap 
functions in the optical limit; both approaches disagree slightly, with the latter having 
problems estimating the total inelastic cross-sections. For the ratio of multiplicities at 
mid-rapidity between the two energies, ^/s NN = 200 GeV and 130 GeV, the experiments 
report values of 1.14 ±0.05 (PHOBOS & BRAHMS), 1.22 ±0.08 (STAR), and 1.17±0.03 
(PHENIX) with no indication of any significant centrality dependence || |1],|],|5||J . 
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Figure 1 : Charged particles per participant pair as 
a function of number of participants for ^/s NN = 130 
and 200 GeV for all four RHIC experiments (from 
|l]|§). PHENIX and STAR data are preliminary. 
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Figure 2: Charged particles per participant pair 
as a function of number of participants in ^/s NN = 
19.6, 130 and 200 GeV Au+Au collisions measured 
by PHOBOS (from §). The curves are two- 
component fits described in the text and in [Q. 



The interpretation of the scaling of the 
multiplicity at mid-rapidity function 
of N part appears still ambiguous. A sim- 
ple model by Kharzeev and Nardi (KN) ex- 
plains the dependence in a two-component 
approach differentiating between soft pro- 
cesses scaling with N part and hard pro- 
cesses scaling with Nb in @. When fit to 
the data as shown in Fig. |2] by the PHO- 
BOS collaboration the model allows one to 
extract the fraction of particles produced 
from hard processes and one obtains val- 
ues of 36% for v /s^ = 130 GeV, and 45% 
for 200 GeV, respectively [f|J§. A second 



class of calculations is based on parton sat- 
uration ||[nj; since the parton densities 
in the initial stage of the collision can be 
related to the density in the final state a 
parametrized dependence of the saturation scale Q s on ^ and impact parameter allows 
one to predict dN^/dr]. However, the predictions from these and related models have 
been found to be almost indistinguishable when applied to RHIC data, especially because 
of the large experimental uncertainties in the calculation of N part for very peripheral colli- 
sions, the region where the differences between the various models become more apparent. 
An exception are models based on final state saturation who significantly overpredict the 
yield at low N part . The systematic uncertainties in N part for peripheral events can only 
be reduced when data from collisions of light ions, A< 100, become available. 
A very surprising finding in the context 
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of particle production was presented by the 
PHOBOS collaboration at this conference 
|||nj. The authors noted that the total 
charged multiplicity per participant in ele- 
mentary e + e~ and central AA collisions is 
identical over a wide range of \fs. That 
alone is a very remarkable fact but they 
also found that the same holds for the pp 
(pp) collisions when compared to the ef- 
fective center-of-mass energy y/s&, that is 
the nominal energy minus the energy car- 
ried away by the leading protons. The au- 
thors actually used ^/s^ = yfs/2 which 
was verified by PYTHIA simulations. This 
is depicted in Fig. [| Also shown is a per- 

turbative QCD calculation for the multiplicity in e + e~ fit to the data. 

The agreement between the three fundamentally different collision systems is remark 
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Figure 3: Total charged multiplicity per participant 
pair as a function of y/s for e + e~, pp (pp), and AA 
collisions. The pp (pp) data is plotted versus y^eff = 
(from p|). 
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200 GeV. One might speculate if this agreement is plainly accidental or if it possibly 
points to a kind of universality in the production of multihadronic final states in high 
energy collisions of any elementary particles over some range in y^i. It is even more sur- 
prising that this should also apply for heavy ion collisions where we have indications that 
the system thermalizes and evolves on a time scale significantly larger than in elementary 
collisions and where the majority of particles are formed late. In comparison, the hadron 
production in e + e~ is mostly from hard gluon radiation (2-,3-and 4-jet events), while in 
pp one observes a mixture of soft and hard processes. 

To shed more light on this apparent sim- 
ilarity it is instructive to look at the scal- 
ing of the average transverse momentum of 
the produced particles in all three systems 
with y/s Any universality in the pro- 
duction mechanism should, to some degree, 
also manifest itself in the evolution of (pt)- 
However, this seems not to hold as illus- 
trated in Fig. | ||. The figure shows that 
the yfs dependence of (p?) in e + e~ is con- 
siderably steeper than in pp (pp), mainly 
due to the absence of soft processes in e + e~ 
collisions. This suggests that the agree- 
ment in particle multiplicity might be ac- 
cidental. Still, the fact remains that the 
total multiplicity per participant appears 
to be remarkably similar for e + e~, nucleon- 
nucleon, and nucleus-nucleus collisions and 
further studies might help to gain insight 
into multihadron production in high energy 
collisions. 

Another interesting finding depicted in 
Fig. [| is the relatively small increase of (py) 
between ^/s NN = 130 and 200 GeV of only 
~ 1%, as pointed out by the STAR collaboration 

hydrodynamics predict a considerably stronger dependence, usually (prY oc dN ch /dr). 
The solid curve is a prediction based on the saturation model, constraint by pp (pp) 
results [R|. The 200 GeV data point clearly falls below this prediction suggesting a 
flattening of the (pt) energy dependence at RHIC. The difference is significant since the 
systematic errors for the 130 and 200 GeV data points are correlated. However, in order 
to prove or falsify any model it is essential to perform further measurements at energies 
between 20 and 200 GeV to study the scaling of (pt) in greater detail. 

Another important observable for characterizing the global properties of bulk matter 
is the transverse energy Et- This was studied in detail by the PHENIX collaboration 
for y/s NN = 130 and 200 GeV [|TJ in the mid-rapidity region. They find that dE T /di] and 
dN c h/dr] increase with N part in a very similar fashion resulting in an almost constant ratio 
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Figure 4: Mean transverse momentum versus y/s 
for e + e~, pp (pp), and AA collisions (from fl6|,p"2|). 
In e + e~ the (pr) is measured with respect to the 
thrust axis. The curve through the e + e~ data points 
are JETSET predictions, the curve through the AA 
data (NA49 and STAR at 130 GeV) is from CGC pre- 
dictions (see text), and that through the pp (pp) data 
is a phenomenological parameterization in log(y / s). 
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(E T )/(N ch )~ 0.9 GeV. This holds for 



130 and 200 GeV. Even more surprising 
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is the fact that studies from Au+Au collisions at y/s NN = 4.8 and Pb+Pb collisions at 
17.2 GeV yield very similar values, suggesting that the increased energy put into the 
system results solely in an increased particle production leaving the average energy per 
particle almost constant. This is depicted in Fig. [5| From the measured dEx/drj for the 
2% most central Au+Au collisions at 200 GeV the authors estimated the Bjorken energy 
density to be £bj ~ 5.5 GeV/fm 3 , assuming a conservative formation time of r = 1 fm/c. 
Similar studies at SPS in Pb+Pb collisions at -Js NN = 17.2 GeV give Ebj ~ 3.2 GeV/fm 3 



14] . These values represent of course only a lower limit for the initial energy density 
since the longitudinal expansion of the system reduces the transverse energy considerably. 
Recent lattice results on QCD thermodynamics estimate the criticial energy density to 
be e ~ 0.70 ± 0.35 GeV/fm 3 [TJ|], a value significantly surpassed already at SPS. 

It is interesting to compare the scaling 
of (3r) I (A?" c h) with that of (pr) depicted in 
Fig. I]. While the small increase in (pt) be- 



tween 



130 and 200 GeV is reflected 
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by the corresponding (Et) / (N^) remain- 
ing constant, the considerably lower value 
of (pt) at CERN/SPS energies appears to 
be contradicted by the transverse energy 
per particle. This, however, could be ac- 
counted for by the different particle com- 
position at SPS and RHIC, although it still 
needs to be verified in quantitative stud- 
ies. It is intriguing to compare the uni- 
versality in (Et) I (iV c h) with an observa- 
tion presented in a paper by Cleymens and 
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Figure 5: dE T /di]\ n=0 /dN ch /dr)\ v= o versus ^/s^ 
for 5% most central events at AGS, SPS, and RHIC 
(from §). 



Redlich in 1998 [|16j in which the authors 
show that the chemical freeze-out param- 
eters obtained at SPS, AGS, and SIS all 

correspond to a unique value of ~ 1 GeV for the average energy per hadron in the local 
rest frame of the system independent of beam energy and mass number. From what we 
have learned so far this empirical observation holds still at RHIC energies leading to a 
considerable unification in the description of hadronic final states in high energy nuclear 
collisions. 



3. Particle Ratios and Chemical Freeze-Out 

One of the most important issues in the physics of heavy ion collisions is the question 
if, and if so at what stage, the produced system thermalizes and to what extent a thermal 
description is appropriate for the evolving system. In order to discuss an equation of 
state and a true order to any associated phase transitions, we need to describe the system 
in terms of a few thermodynamic properties. The use of thermodynamic concepts to 
describe multi-particle production has a long history beginning with Hagedorn in the 
early 1960's |17[]. The concept of a temperature applies, strictly speaking, only to systems 
in at least local thermal equilibrium. Thermalization is normally only thought to occur 
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Figure 6. Particle ratios for 200 GeV central Au+Au collisions at RHIC |,|8|JT|,|^,|l]j22|j2| . The error 
bars depict the statistical errors, the brackets the systematic uncertainties. In the upper row the vertical 
solid lines indicate the RHIC average, the dashed lines the error on the mean. Open circles are used for 
data points which were extracted from figures. Results are preliminary, except BRAHMS and PHOBOS 
data in upper row. 



in the transverse degrees of freedom as reflected in the Lorentz invariant distributions of 
the particles. The measured hadron spectra contain two pieces of information: (i) their 
normalization, i.e. their yields and ratios, provide the chemical composition of the fireball 
at the chemical freeze-out point and (ii) their transverse momentum spectra which provide 
information about thermalization of the momentum distributions and collective flow. It 
is obvious that the observed single particle spectra do not reflect earlier conditions, i.e. 
the hot and dense deconfined phase, where chemical and thermal equilibrium may have 
been established, since rescattering erases most traces from the dense phase. Only those 
effects which are accumulative during the expansion, such as flow, remain. 

The assumption of a locally thermalized source in chemical equilibrium can be tested by 
using statistical thermal models to describe the ratios of various emitted particles. This 
yields a baryon chemical potential /xg, a strangeness saturation factor 7 S , and the tem- 
perature T c h at chemical freeze-out. Because of the absence of any dynamic assumptions 
many details can never be fully absorbed by these models. Discrepancies between model 
and data up to 30% should be considered inside the systematic uncertainty of the thermal 
model approach [p4 |. So far these models are remarkably successful in describing parti- 
cle ratios at SPS f25,26f and now also at RHIC |27|13. This observat ion, together with 
the large collective flow (radial and elliptic) measured at RHIC, is generally considered a 
strong hint that chemical equilibrium is indeed reached. The wide reaching implications 
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Figure 7. Left panel: comparison between RHIC experimental particle ratios for ^/s NN = 130 GeV 
and statistical model calculations with T c h = 176 MeV and jiB — 41 MeV (from f p7pl| |). Right panel: 
comparison between RHIC ratios at v /s NN = 200 GeV and prediction discussed in the text (also f^Jjlf ) . 



of thermal statistical models and the models themselves were the subject of a dedicated 
podium discussion at the conference. I'm not going to summarize the issues brought up in 
this discussion but rather will concentrate on the recent experimental results on particle 
ratios and their implications in the framework of the models. 

Figure ^| shows a compilation of most results on particle ratios at mid-rapidity for 
y/s NN = 200 GeV presented at the conference [^, p^ , p!9| , p0| , ^l] , p2] , ^| . Given the good agree- 
ment among the experiments for the identical particle ratios and the high level of quality 
of the data, it is tempting to calculate the RHIC averages for these results. Adding the 
statistical errors and the systematical uncertainties in quadrature to derive the weights 
one obtains: tt-/tt + = 1.02 ± 0.02, K~/K + = 0.95 ± 0.03, and p/p = 0.75 ± 0.03. All 
experiments reported to observe no pt or centrality dependence of these ratios for pt < 3 
GeV/c, confirming earlier results at 130 GeV. STAR reports a decrease in p/p at high p?; 
for discussion on this topic and also on the significant increase of the p/ir ratios from low 
to medium p T see |2^J3^ . 

The analyses on non-identical particle ratios are not as complete and in most cases 
the systematic uncertainties are still under evaluation. While the net-baryon chemical 
potential at chemical freeze-out is essentially determined by the baryon to antibaryon 
ratios (p/p, A/A etc.), the non-identical particle ratios are the "thermometer" of the ther- 
mal statistical models. Because of the lack of sufficient constraints it would be therefore 
premature to invoke the thermal model on the 200 GeV data at this point. It is, how- 
ever, instructive to compare the currently available ratios with predictions made in p7 |. 
Here the authors used a phenomenological parametrization of yUe, obtained from thermal 
parameters derived from the statistical model at lower energies in conjunction with the 
assumption of constant energy per particle (see above) to extrapolate to 200 GeV. 

This comparison, together with the statistical model fit for 130 GeV, updated with the 
latest values presented at this conference is shown in Fig. |7| ||31|| . The predictions match 
well with the current results and indicate no significant change in T c h but a drop in \xb 
from ~ 41 MeV at v /s^7 = 130 GeV to 29 MeV at 200 GeV. The latter value is also 
in agreement with calculations made by various authors at the conference using different 
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approaches H5| Jl9|Jl8| . The chemical freeze-out temperature is naturally limited by the 
confinement phase-transition temperature assumed to be around 175 MeV, although T c h 
is actually not constrained in thermal model fits. 

Another interesting result was reported 
by the BRAHMS collaboration who pre- 
sented a detailed study of identical particle 1 — •• J- 

ratios function of rapidity for central 
events |ggO 
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As shown in Fig. | 
ratio is consistent with unity over the con- 
sidered rapidity range while the K~/K + 
ratios drops by ~ 30% at y = 3 from its 
mid-rapidity value and the p/p ratio by 
~ 70%. Interestingly, all ratios remain con- 
stant for \y\ < 1, consistent with the as- 
sumption of boost invariance around mid- 
rapidity (see below). From those data the 
authors derived a net-baryon chemical po- 
tential at y = 3 of \ib ~ 130 MeV within 
the framework of a statistical model, assuming that the particle sources in the different y 
regions are in local chemical equilibrium and that strangeness is locally conserved. How- 
ever, to what extent a thermal interpretation at large forward rapidities is justified is 
subject to further studies. 



Figure 8: Antiparticle to particle ratios as a func- 
tion of rapidity in 200 GeV Au+Au collisions (from 



4. Boost Invariance at RHIC? 



Many models assume directly or indi- 
rectly that the system created in heavy 
ion collisions is boost invariant, i.e. invari- 
ant under Lorentz transformations in the 
beam direction. Commonly, but not cor- 
rectly, one considers a system boost in- 
variant within a given rapidity interval if 
the rapidity distribution dN/dy is constant 
within that range. Strictly speaking, how- 
ever, it requires all Lorentz invariant ob- 
servables to remain constant. The pseu- 
dorapidity distribution of charged particles 
at RHIC shows a plateau extending over 
almost 3 units [[32]. dN/dr] distributions, 
however, can not be used to draw any con- 
clusions on boost invariance because the 
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Figure 9: Rapidity distributions for 7r ± , K ± ,p, and 
p in 200 GeV Au+Au collisions for the top 10% most 
central events (data taken from |5jj2(]|). 



Jacobian dy/dr](pT,v) that tends to flatten otherwise peaked rapidity distributions. The 
BRAHMS collaboration presented at this conference the rapidity distribution of pions, 
kaons, and protons and their antiparticles over 6 units of rapidity < 3) for Au+Au 



collisions at 



200 GeV 
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As shown in the compilation of all six 
distributions in Fig. |9] and given the sta- 
tistical errors and systematic uncertainties 
a plateau is at most limited to ±1 unit in 
rapidity. The upper panel in Fig. |l^ shows 
the pion rapidity distribution in this very 
region as measured by STAR with higher 
granularity |J . A study of the slope of the 
corresponding p T -distribution, depicted in 
the lower panel, shows that the slope starts 
to decrease from its mid-rapidity value at 
around ±i] ~ 0.5 thus limiting the region 
of boost invariance to \rj\ < 0.5. This is 
also confirmed for the case of the protons 
by BRAHMS where no significant change 
in slope is observed between between y — 
and 0.9 [20]. In this context it is also im- 
portant to recall Fig. § to verify that the 
particle to antiparticle ratios of the most 
abundant particles (n, K, and p) remain 
constant within half a unit around mid- 
rapidity. Recent studies of elliptic flow also indicate a continous drop in v 2 for \r}\ > 1 that 
is, at least in parts, related to the strong dependence of Vi onpT(rj) f33f . We conclude that 
the currently available RHIC data suggests a rather small truely boost-invariant region 
of at most \r)\ < 0.5. 
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Figure 10: Upper panel: 7r~ rapidity distribution 
in 130 GeV Au+Au collisions for top 5% most cen- 
tral events. Lower panel: inverse slope obtained from 
Bose-Einstein fits to the pt spectra at different ra- 
pidities (from [||,D). 



5. Inclusive Transverse Momentum Spectra and Radial Flow 

Inclusive hadron spectra were intensively discussed at this conference in the context of 
the observed suppression of high-p^ yields for central collisions when compared to either 
peripheral collisions or pp (pp) reference data [j30fl . These studies focus naturally only 
on a tiny fraction of all produced particles. The majority of particles emitted from the 
system are soft. At RHIC 99.9% of all charged particles have momenta below 2 GeV/c, 
far outside the range of perturbative QCD. 

Transverse momentum spectra of identified particles reflect the system at kinetic freeze- 
out and allow us to extract information from the latest stage of the evolution when 
the system was still thermally coupled and governed by elastic interactions among its 
constituents. The measured inverse slope parameter is determined by two components: 
the actual temperature at the freeze-out and the transverse flow component. In simple 
terms this can be approximated as T = Tf Q + m(/?T) 2 where Pt is the transverse flow 
velocity. This ansatz, however, has the disadvantage that attempts to extract Tf Q and 
Pt are strongly dependent on the range in which the slopes were determined. Of even 
greater concern is the assumption of a fixed flow velocity which oversimplifies the problem 
considerably. To overcome these problems and to avoid the complexity of adjusting the 
initial energy density and the equation of state in a full hydrodynamical model calculation 
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many studies now use the so-called 'blastwave' parametrization [34]]. Here, the invariant 
cross-section is fit to 
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where p = tanh -1 P T is the transverse rapidity and (3? = (3 s (r / R max ) n depends on the 
chosen flow profile and the flow at the surface (3 S . There is no commonly accepted flow 
profile and n varies in the different analysis between 0.5 and 2. It is important to keep in 
mind that Tf Q and j3r are correlated. Increasing Tf Q or (3r has to some degree a similar effect 
on the spectral shape. This problem can be overcome by applying a blastwave motivated 
parametrization also to HBT radii, R(kj-), since here Tf Q and (3t are anti-correlated helping 
to further constrain T fo and thus /3y. 
Van Leeuwen presented an im- 
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pressive compilation of NA49 
transverse mass spectra in 40, 
60, and 158 AGeV collisions 
Figure |11| depicts the re- 
sults from Pb+Pb collisions at 
158 AGeV where all spectra were 
fit to a blastwave parameteriza- 
tion. The obtained freeze-out pa- 
rameters are Tf Q = 122 — 127 
MeV and {(3 T ) = 0.48 ± 0.01 in 
agreement with previous studies. 
Blastwave fits to RHIC data give 
a slightly smaller freeze-out tem- 
perature Tf ~ 110 MeV but a 
higher flow value (/3r) = 0.55 — 
0.6 due to the higher pressure in 
the system The surprising 

finding in the NA49 analysis, how- 
ever, is the fact that not only are 

the 7r, K, p, and A spectra well described by the fit but so are the spectra of the multi- 
strange baryons 5 and fl So far multi-strange baryons were assumed to show less or no 
flow, due to a possibly very small elastic cross-section. This was supported by simple fits 
to the Q spectra which yield Tf Q values close to the chemical freeze-out temperature. It 
has been speculated that if the elastic cross- sect ions for the Q were indeed very small, 
a non-zero flow component could be interpreted as a signature for partonic flow. The 
situation at RHIC is still ambiguous since the current Q mr-spectra lack statistics and 
within the errors can be fit by either assuming Tf Q ~ 170 and (/3t) ~ 0.4 or Tf Q = 130 and 
(Pt) ~ 0.5 [57|. One can, however, already exclude the cases of no flow and, unlike at 
SPS, it cannot be described by parameters obtained from fits to pions, kaons, and protons 
alone. To what extent a combined fit of all spectra, similar to that presented for the SPS, 
yields better results still has to be seen. Higher precision data is certainly needed here. 

K~ 



Figure 11: Transverse mass spectra in central 158 AGeV 
Pb+Pb collisions measured by NA49. The curves show the 
result of blastwave fits with a constant flow profile, i.e. n = 
(from @). 



Figure [I2| shows a compilation of pr spectra for n 



and p from all four RHIC 
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p T (GeV/c) 

Figure 12. Compilation of preliminary transverse momentum spectra of tt~,K~, and p for 200 GeV 
central Au+Au collisions from all RHIC experiments. Except for the PHOBOS data all spectra are not 
feed-down corrected. The curves are fits to PHOBOS, BRAHMS, and PHENIX data. For n' a power 
law in rriT was used; Maxwell-Boltzmann distributions for K~ and p. 

experiments at ^/s^ = 200 GeV. For the first time the PHOBOS collaboration presented 
data at very low pt, down to 30 MeV for pions, thus extending our knowledge into a range 
which is not accessible to the other three experiments. The low-pr region is very sensitive 
to dynamic effects and will help to constraint models even further. Except for PHOBOS 
all spectra shown in the figure are not feed-down corrected. It is important to note, that 
the RHIC experiments are not equally sensitive to feed-down, e.g. BRAHMS because of its 
small aperture is less affected than a large acceptance detector such as STAR. Note, also 
the slightly different centrality selections. The difference in yields between 15% and 5% 
centrality is approximately 15%. One of the problems in determining the absolute yield 
of particles is the need to extrapolate to pt = 0. At RHIC this results in considerable 
uncertainties in the yields of typically ~ 5 — 10% for experiments limited to pt > 200 
MeV/c or above. The new low-pr data now allows us to significantly reduce this source of 
uncertainty. It is interesting to note that pion spectra actually do not follow a power-law 
in pt down to low px as indicated by earlier measurements, but can be very well described 
by a power-law in ttit'- oc A (1 + m^/mo) - ™ as indicated by the fit shown in Fig. |P2| . 

Radial flow affects the spectra of heavier particles considerably more than light particles. 
One of the consequences of strong radial flow at RHIC is that the yields of pions, kaons, 
and protons essentially become equal at around pt = 2 GeV/c as depicted in the smaller 
panel in Fig. [12]. At large pt the baryon to meson ratios should drop again as flow should 
affect high-pr particles to a much lesser degree and the ratios should approach values 
predicted by pQCD. Other studies follow a different approach, invoking novel baryon 
dynamics attributed to gluonic baryon junctions that predict the baryon-enhancement 
only in a finite moderate-pr window ||38|| . However, the 'turn-over' point, predicted in 
these models to be around 3 GeV/c, has not yet been clearly identified [|T^] . 



12 



6. Summary 

The recent results highlighted at the conference and in this review show the enormous 
efforts of the community to study 'bulk' matter governed by the non-perturbative regime 
of QCD. The new RHIC data recorded at y/s NN = 200 GeV seem, without exceptions, 
to confirm the picture that evolved from the first data, at 130 GeV, reported at Quark 
Matter 2001. This picture cannot be understood in terms of global variables, ratios, and 
spectra alone but only by combining the information we get from studies of HBT, elliptic 
flow, and high-p^, to name only a few. 

Although we observe quantitative differences between SPS and RHIC in many parame- 
ters, there appears not to be any striking qualitative difference between these two energy 
regimes with the very prominent exception of the onset of hard scattering processes at 
RHIC. 

The study and interpretation of soft hadron spectra, particle ratios, and yields is domi- 
nated by the success of the thermal statistical models. Although critical issues concerning 
these models need to be resolved, it remains a fact that they describe the ratios and yields 
remarkably well over a wide range of energies. They indicate that in the energy range from 
a few GeV up to RHIC energies the observed hadrons originate from a system in chemical 
equilibrium along a unified freeze-out curve. This curve provides the relation between the 
temperature and the baryon chemical potential. At RHIC we are approaching an almost 
net-baryon free system with fig ~ 25 MeV where, for the first time, more baryons are 
produced at mid-rapidity than transported from beam-rapidity. The chemical freeze-out 
temperatures at SPS (T c h ~ 165 MeV) and at RHIC (T c h ~ 175 MeV) extracted from 
these models appear to be close to the critical temperature. This implies that chemical 
equilibrium is not caused by kinetic equilibration through hadronic rescattering but in- 
dicates that hadron formation proceeds by statistical hadronization from a prehadronic 
state. 

Detailed analysis of the hadron spectra shows that the system expands collectively 
under strong internal pressure. Radial flow at RHIC appears to be slightly higher than at 
SPS, the kinetic freeze-out temperatures are very close to each other, possibly somewhat 
lower at RHIC (T fo ~ 110 MeV) than at SPS (T fo ~ 120-130 MeV). Studies of resonances 
[ [39(1 and correlations |^(| show that the time scale for emission is very short (2-3 fm/c) 
while the overall lifetime of the system appears to be on the order of 10 fm/c. This implies 
that after hadronization the hadron abundances freeze out more or less immediately. 

We could declare proof of the quark-gluon plasma on the basis of indirect evidence, but 
the fact that a new phase exists is almost trivial compared to characterizing its features. 
The interpretation of bulk properties in heavy ion collisions was, and still is, complex. 
We have evidence for thermalization of hot and dense matter and we have indications of 
unusual behavior in rare, high momentum probes. The level of collectivity is surprising 
but the timescales are puzzling. We observe matter that is surely not a simple collection 
of elementary particles, and we have the tools to study it. 

Let me finally thank all the speakers at the conference who have provided me with 
their results. I am grateful to many colleagues for critical discussions and useful support, 
in particular to M. Baker, R. Bellwied, M. Calderon, H. Caines, J. Harris, B. Hippolyte, 
F. Laue, D. Magestro, K. Redlich, C. Suire, and Z. Xu. 
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